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Fractional pure birth processes 
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We consider a fractional version of the classical nonlinear birth process of which the Yule-Furry 
model is a particular case. Fractionality is obtained by replacing the first order time derivative 
in the difference-differential equations which govern the probability law of the process with 
the Dzherbashyan-Caputo fractional derivative. We derive the probability distribution of the 
number M v if) of individuals at an arbitrary time t. We also present an interesting representation 
for the number of individuals at time t, in the form of the subordination relation N v (t) = 
N(T2 V (t)), where N{t) is the classical generalized birth process and T2 V {t) is a random time 
whose distribution is related to the fractional diffusion equation. The fractional linear birth 
process is examined in detail in Section 3 and various forms of its distribution are given and 
discussed. 

Keywords: Airy functions; branching processes; Dzherbashyan-Caputo fractional derivative; 
iterated Brownian motion; Mittag-Leffler functions; nonlinear birth process; stable processes; 
Vandermonde determinants; Yule-Furry process 

1. Introduction 

We consider a birth process and denote by N{t), t > 0, the number of components in 
a stochastically developing population at time t. Possible examples are the number of 
particles produced in a radioactive disintegration and the number of particles in a cosmic 
ray shower where death is not permitted. The probabilities Pk(t) = Pr{AT(t) = k} satisfy 
the difference-differential equations 

= -\ k Pk + A fc _ip fc _i, k>l, (1.1) 

at 

where, at time t = 0, 

pm={1 III (i.2) 

This means that we initially have one progenitor igniting the branching process. For 
information on this process, consult Gikhman and Skorokhod [5], page 322. 



This is an electronic reprint of the original article published by the 1SI/BS in Bernoulli, 
2010, Vol. 16, No. 3, 858-881. This reprint differs from the original in pagination and 
typographic detail. 



1350-7265 © 2010 ISI/BS 



Fractional birth processes 



859 



Here, we will examine a fractional version of the birth process where the probabilities 
are governed by 

= ~\ k Pk + \k-iPk-i , k>l, (1.3) 

at v 

and where the fractional derivative is understood in the Dzherbashyan-Caputo sense, 
that is, as 

d»p k 1 f* (d/ds)p fc («) , 



for < v < 1 (1.4) 

dt" r(i - v) y (t - s) v y ' 

(see Podlubny [12]). The use of a Dzherbashyan-Caputo derivative is preferred because 
in this case, initial conditions can be expressed in terms of integer-order derivatives. 

Extensions of continuous-time point processes like the homogeneous Poisson process 
to the fractional case have been considered in Jumarie [7], Cahoy [3], Laskin [9], Wang 
and Wen [17], Wang, Wen and Zhang [18], Wang, Zhang and Fan [19], Uchaikin and 
Sibatov [15], Rcpin and Saichev [13] and Bcghin and Orsingher [2]. A recently published 
paper (Uchaikin, Cahoy and Sibatov [16]) considers a fractional version of the Yule-Furry 
process where the mean value EN v (t) is analyzed. 

By recursively solving equation (1.3) (we write Pk(t), t > 0, in equations (1.3) and p^{t) 
for the solutions), we obtain that 

(1.5) 

( fe-l k ( n 

II A ; E TT k — —E v>1 (-X m t v ) , k > 1, 

j=l m=l^LLl=l,l^m\ A i Am > ) 

{E vA {-X 1 t u ), k = l. 

Result (1.5) generalizes the classical distribution of the birth process (see Gikhman and 
Skorokhod [5], page 322, or Bartlett [1], page 59), where, instead of the exponentials, we 
have the Mittag-LefHcr functions, defined as 

oo h 

g MW = E r J + iv xeK,^>0. (1.6) 

h=0 ^ ' 

The fractional pure birth process has some specific features entailed by the fractional 
derivative appearing in (1.4), which is a non-local operator. The process governed by 
fractional equations (and therefore the related probabilities p v k it) = Pr{N v (t) = k}, k > 1) 
displays a slowly decreasing memory which seems a characteristic feature of all real sys- 
tems (for example, the hcrcditariety and the related aspects observed in phenomena such 
as metal fatigue, magnetic hysteresis and others). Fractional equations of various types 
have proven to be useful in representing different phenomena in optics (light propagation 
through random media), transport of charge carriers and also in economics (a survey of 
applications can be found in Podlubny [12]). Below, we show that for the linear birth 
process N v (t),t>0, the mean values EiV„(t), V&iN v (t) are increasing functions as the 
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order of fractionality v decreases. This shows that the fractional birth process is capa- 
ble of representing explosively developing epidemics, accelerated cosmic showers and, in 
general, very rapidly expanding populations. This is a feature which the fractional pure 
birth process shares with its Poisson fractional counterpart whose practical applications 
have been studied in recent works (see, for example, Laskin [9] and Cahoy [3]). 
We are able to show that the fractional birth process M v (t) can be represented as 

N v {t)=N{T 2v {t)), t>0,0<i/<l, (1.7) 

where T% v (t), i > 0, is the random time process whose distribution at time t is obtained 
from the fundamental solution to the fractional diffusion equation (the fractional deriva- 
tive is defined in (1.4)) 

d 2v u d 2 u n .„ . 

subject to the initial conditions u(s, 0) = S(s) for < v < 1 and also ut(s,0) = for 
1/2 < v < 1, as 

¥r{T 2v {t)eds} = { 2u ^ s ^ As fOTS >° ; (1.9) 
L for s < 0. 

This means that the fractional birth process is a classical birth process with a random 
time T2„(i) which is the sole component of (1.7) affected by the fractional derivative. 
In equation (1.8) and throughout the whole paper, the fractional derivative must be 
understood in the Dzherbashyan-Caputo sense (1.3). The representation (1.7) leads to 

/>oo 

Pi{/V v (t) = k} = / Pi{iV(s) = k}PT{T 2 „(t)eds}, (1.10) 

JQ 

where 

H Aj Snt Wm (A i -A„ l )' fc>1 ' s>0 ' (i.ii) 

e" AlS , fc = l, s>0. 

Formula (1.10) immediately shows that J2k P r {A/^(t) = k} = 1 if and only if 
^ fe Pr{A/"(t) = k} = 1. It is well known that the process N(i), t > 0, is such that 
Pr(Af(t) < oo) = 1 for all t > (non-exploding) if ^2 k X^ 1 = oo (see Feller [4], page 
452). 

A special case of the above fractional birth process is the fractional linear birth process 
where = Xk. In this case, the distribution (1.5) reduces to the simple form 

pm =E ( * Z I) i-iy-'EvA-vn k > i,t > o. (i.i2) 



Fractional birth processes 



861 



For v = 1, we retrieve from (1.12) the classical geometric structure of the linear birth 
process with a single progenitor, that is, 

^(t) = (l-c- At ) fc - 1 c- At , fc>l,t>0. (1.13) 

An interesting qualitative feature of the fractional linear birth process can be extracted 
from (1.12); it permits us to highlight the dependence of the branching speed on the 
order of fractionality v. We show in Section 3 that 

Pr{A^(di) = n + 1|JV„(0) = no} ~ ( L14 ) 

r(w + i) 

and this proves that a decrease in the order of fractionality v speeds up the reproduction 
of individuals. We are not able to generalize (1.14) to the case 

Pr{iV„(t + dt) = n + l|JV„(t) = n } (1.15) 

because the process we are investigating is not time-homogeneous. For the fractional 
linear birth process, the representation (1.7) reduces to the form 

N v (t) = N{T 2v (t)), f>0,0<i/<l, (1.16) 

and has an interesting special structure when v = 1/2". For example, for n = 2, the 
random time appearing in (1.16) becomes a folded iterated Brownian motion. This means 
that 

N 1/i (t) = N(\B 1 (\B 2 (t)\)\). (1.17) 

Clearly, |2?2(*)| is a reflecting Brownian motion starting from zero and |£>i(|£>2(i)|)| is a 
reflecting iterated Brownian motion. This permits us to write the distribution of (1.17) 
in the following form: 



Pr{JV 1/4 (*) = fe|iV 1/4 (0) = 1} 

> C roc -s 2 /(4w) -tu 2 /(4t) 

I Jo V%n2u V2n2t 



(1.18) 



The case v = 1/2™ involves the (n — l)-times iterated Brownian motion 

l n . 1 (t)=B 1 (\B 2 (---\B n (t)\ ■••)!) (1.19) 

with distribution 

Pr{|B 1 (|B 2 (...|B n (*)|...))|€d S } 



= ds2™ / — -j^=^^- dwi / — -j^=^^- Auji ••• I : — du ri 

Jo v4ttwi Jo v47tw 2 J V47rf 

For details on this point, see Orsingher and Beghin [11]. 



(1.20) 
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2. The distribution function for the generalized 
fractional birth process 

We now present the explicit distribution 

Pr{Af v (t) = k\Af v (0) = l}=pl{t), t>0,k>l,0<v<l, (2.1) 

of the number of individuals in the population expanding according to (1.3). Our tech- 
nique is based on successive applications of the Laplace transform. Our first result is the 
following theorem. 

Theorem 2.1. The solution to the fractional equations 



AfcPfe + \k-1Pk-1, k>l,0<v<l, 



1, k = l, 
0, fc>2, 



Pfc(0) 



is given by 

p%(t) = -Pi{JV„(t) = k\JV 1/ (0) = l} 

1 




In view of the fact that 

H v + X 



(2.2) 



(2.3) 

—-E^(-x m n}, k>i, 

Eu^-XiV), k = l. 

Proof. We prove the result (2.3) by a recursive procedure. 
For k = 1 , the equation 

& = -A lPl) Pi(0) = l, (2.4) 

is immediately solved by 

p v 1 {t) = E v , l (-\ l t v ). (2.5) 
For k = 2, equation (1.3) becomes 

P2(0)=0. 



e-^ I i(-Ait")di=-J- ri (2.7) 
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the Laplace transform of (2.6) yields 

r ( \ 



1 



1 



[i v + Ai ^ + A 2 



A 2 — Ai 

In the light of (2.7), from (2.8), we can determine the probability p\{£): 

Ai 



P » 2 (t) = [E Uil (-\ 1 t»)-E vA (-\ 2 t v )] 



A2 — Ai 



Now, the Laplace transform of 



A3P3 + r^-^iC-Aif) - ^.iC-Aaf)] 



d^ ^ A 2 - Ai 

yields, after some computation, 

L 3 (/z) = X 2 Xin" 



1 



1 



(A 2 -A 1 )(A 3 -Ai)^ + Ai 

1 1 

+ (Ai-A 2 )(A 3 -A 2 )^ + A 2 

1 1 

+ (Ai-A 3 )(A 2 -A 3 )^ + A 3 



From this result, it is clear that 
p v 3 (t) = X 2 X 1 



1 



(A 2 -A 1 )(A 3 -A 1 ) 
1 

+ (Ai-A 2 )(A 3 -A 2 ) 
1 



(Ai-A 3 )(A 2 -A 3 ) 



E v>1 (-\ 3 t») 



The procedure for k > 3 becomes more complicated. However, the special 
instructive and so we treat it first. 

The Laplace transform of the equation 



d"p 4 



= — XiPi + AiA 2 Ag 



1 



(A 2 -Ai)(A 3 -Ax) 
1 

+ (Ai-A 2 )(A 3 -A 2 ) 
1 



(Ai-A 3 )(A 2 -A 3 ) 



^,i(-Aif) 

E^(-\ 2 n 

E vA {-X 2 t v ) 
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subject to the initial condition p 4 (0) = 0, becomes 

1 



L 4 (/i) = AiA 2 A 3 // 1 ' 



(A 2 - Ai)(A 3 - Ai)(A4 - Ai) {/i" + Ai fi v + A 4 



1 



(Aj - A 2 )(A 3 - A 2 )(A 4 - A 2 ) \ ii v + A 2 //" + A 4 



1 



1 



1 



1 



1 



(Ai - A 3 )(A 2 - A 3 )(A 4 - A 3 ) i fi v + A 3 fi v + A 4 
The critical point of the proof is to show that 

-((A 3 - A 2 )(A 4 - A 2 )(A 4 - A 3 ) - (A 3 - A1XA4 - Ax)(A 4 - A 3 ) 
+ (A 2 -A 1 )(A 4 -A 1 )(A 4 -A 2 )) 

1 

X (A 2 - Ai)(A 3 - Ax)(A 4 - A!)(A 3 - A 2 )(A 4 - A 2 )(A 4 - A 3 ) 
1 



(Ai-A 4 )(A 2 -A 4 )(A 3 -A 4 )' 



We note that 



= det 



(1 1 1 1\ 
1111 

Ai A 2 A3 A 4 
\\\ \\ A§ A 4 / 

111 

det I A 2 A 3 A 4 

^2 ^3 , 
111 



111 



111 



det Ai A 3 A 4 + det Ai A 2 A 4 - det A x A 2 A 3 



\ 2 \ 2 x 2 

> A l A 3 A 4 



X 1 A 2 A 4 



, A i ^2 ^3 



= (A 3 - A 2 )(A 4 - A 2 )(A 4 - A 3 ) - (A 3 - Ax)(A 4 - A X )(A 4 - A 3 ) 
+ (A 2 - Ai)(A 4 - Ai)(A 4 - A 2 ) - (A a - Xi)(X s - Ax)(A 3 - A 2 ), 

where, in the last step, the Vandcrmonde formula is applied. 
By inserting (2.16) into (2.14), we now have that 



i 4 (/i) = \1\2\3n" 



1 



1 



(A 2 -A 1 )(A 3 -A 1 )(A 4 -A 1 )^ + A 1 

1 1 

+ (Ai-A 2 )(A 3 -A 2 )(A 4 -A 2 ) M - + A 2 



(2.14) 



(2.15) 



(2.16) 
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(Ai - A 3 )(A 2 - A 3 )(A 4 - A 3 ) fi v + A 3 

1 1 
(Ai - A 4 )(A 2 - A 4 )(A 3 - A 4 ) \i v + A 4 
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so that by inverting (2.17), we obtain the following result: 



4 



j = l (tn 



=i ri/=i,;^m(A; - A m ) 



E v .\ (— \ m t v ) 



(2.18) 



We now tackle the problem of showing that (2.3) solves the Cauchy problem (2.2) for all 
k > 1 , by induction. This means that we must solve 



1 



fc-l ( k-l 
j=l (m=l lll=l,l^m\ A l A ™> 



U*(o) = o, 

The Laplace transform of (2.19) reads 



fc>4. (2.19) 



fe-i 



fe-i 

E 



1 U=l rii=l,(^ m (Ai - Am) M" + Am 



A^- 1 ^ 1 



E 



^+ A ^ 1 nL« m (A i -Am) 



(2.20) 



We must now prove that 



fc-i 
E r 



i Ili=i,z^m(Ai - A m ) Ui=i,i^k( x i - Afc) 

and this relation is also important for the proof of (1.11). 
In order to prove (2.21), we rewrite the left-hand side as 



(2.21) 



y> nLint>h(A( - A h ) 



i nLij#™(Ai - A m ) yiUi nf>ft(Ai - Ah) 



(2.22) 
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and concentrate our attention on the numerator of (2.22). By analogy with the calcula- 
tions in (2.16), we have that 



= det 



/ 1 
1 

Ai 



i 
i 

A-2 
v fc-2 



1 
1 

A?r 2 



Afc 



Ar 2 



i 



i 



i 



^(-lr-Mcti 



m=l 



A TO — i A TO _|_i 



Vn+l 




(2.23) 



Hh=lHl>h( X l - X h) 



fe-1 T-r/c-l 



rutins -A/0 , ntinf> fc (Ai-A A ) 



E 



^i rTi=i,/^ m (A; - A m ) m=1 rii=i,/^ m (Ai - A TO ) ni=i,;^fe(A; - Afc) 

In the third step of (2.23), we applied the Vandcrmonde formula and considered the fact 
that the nth column is missing. It must also be taken into account that 

nf>i(A;-Ai) n/>2(A/-A 2 ) Ilf>m-i( A i _ A m _i) 



(A m — Ai) 



(A m — A2) 



(A m — A m _i) 



Il/>m(A/ - A,: 



T\l>,n( X l ~ A m ) 

l>m+l 

From (2.22) and (2.23), we have that 
fc-i 



' (Aj-A m+1 )--- " [ (A ; -A fc _! 



(2.24) 



;>fc-i 



E 



1 



=1 Ili=i.i^m(A; - A, ; 



g ntlntjA^A,) ^1 

m=l Ili=l,Z^m(Az - Am) Y\ h =l Ul>h _ Afc) 



(2.25) 



]l/=i,i^fc(Ai - A fe ; 
In view of (2.25), we can write that 



fc-i 



ifc(M) = IlA,E 



j = l ra= 



=inLi«™(A ; -A m ) ^ + A, 



(2.26) 



because the fcth term of (2.26) coincides with the last term of (2.20) and therefore, by 
inversion of the Laplace transform, we get (2.3). □ 
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Remark 2.1. We now prove that for the generalized fractional birth process, the rep- 
resentation 

Af v (t)=Af(T 2l/ (t)), t>0,0<i/<l, (2.27) 

holds. This means that the process under investigation can be viewed as a generalized 
birth process at a random time T2„(i), t > 0, whose distribution is the folded solution to 
the fractional diffusion equation (1.8). 



/>oo 

/ e-"*ft,(u,t)di 
Jo 



by L 2.3) n fy n\ ± n M:^\ +u^(-Aio [e-'-dt 

J0 U=2 3 '=1 m=l llj^mW — A »>V J 

= Vw fe TTA-V-^ 1 | ^ (2 28) 

/■oo f oo fc— 1 fc u-\ ~\ 

= E « fc II X i E "HT^ e -^"+A m ) + U e- S ^"+ Al ) d,s 

"' U=2 j=l m=l Ilj^m(Aj - Am) J 

/* OO /*OG />00 

= / 0(u,sy- 1 e- s ""'ds = / £?(«,«) / e-^fejs.tJdMs 
Jo Jo Jo 

e-"*( /°°gK s )/ T2 „( s ,t)d s ldt, 



where 



/•OO 

/ e^/r^MM^^c"^, s>0, (2.29) 
Jo 

is the Laplace transform of the folded solution to (1.8). From (2.28), we infer that 



g v {u,t)= / g(u,s)f T2 As,t)ds (2.30) 
Jo 

and from this, the representation (2.27) follows. 

Remark 2.2. The relation (2.27) permits us to conclude that the functions (2.3) are 
non-negative because 



/>oo 

Pr{M,(*) =k}= Pr{Af(s) = fc} Pr{T a „(*) 6 ds}, 
Jo 



(2.31) 



and Pr{7V(s) = fc} > and ^ fe Pr{7V(s) = fc} = 1, as shown, for example, in Feller 
[4], page 452. Furthermore, the fractional birth process is non-exploding if and only 
if X]fc(VAfc) = o° f° r ah values of < v < 1. 
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3. The fractional linear birth process 

In this section, we examine in detail a special case of the previous fractional birth process, 
namely the fractional linear birth process which generalizes the classical Yule-Furry 
model. The birth rates in this case have the form 

A fc = Afc, A>0, k>l, (3.1) 

and indicate that new births occur with a probability proportional to the size of the 
population. We denote by N u (t) the number of individuals in the population expanding 
according to the rates (3.1) and we have that the probabilities 

pl(t)=Pi{N v (t) = k\N v (0) = l}, k>l, (3.2) 

satisfy the difference-differential equations 

^-^ = -Xk Pk + X(k-l)p k - 1 , 0<u<l,k>l, 

(3.3) 



Pfc(O) 



1, k = l, 
0, k>2. 



The distribution (3.2) can be obtained as a particular case of (2.3) or directly, by means 
of a completely different approach, as follows. 

Theorem 3.1. The distribution of the fractional linear birth process with a simple initial 
progenitor has the form 



(3.4) 

E ( ) ~- 1 ) (-ly-'EuA-vn, k > i,o < v < i, 



pl(t)=Pv{N v (t) = k\N v (0) = l} 

k 

.3=1 v 

where E v ^i(x) is the Mittag-Leffler function (1-6). 

Proof. We can prove the result (3.4) by solving equation (3.3) recursively. This means 
that p u k _Y (t) has the form (3.4), so p^{t) maintains the same structure. This is tantamount 
to solving the Cauchy problem 

lp fc (0) = 0, k>l. 
By applying the Laplace transform Lk^(n) = J °° e _M *pfc(i) dt to (3.5), we have that 



'fe-i 



j - i / py + Xj \py + Xk 
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Conveniently, the Laplace transform (3.6) can be written as 



fj,"+X n v + Xk 
(*-!)(* -2) 



(fc-1) 



1 



1 



1 



ix v + 2A fx v + Afc 
1 



+ (k -1)(-1) 



fc-2 



+ 3A /a" + fcA 
1 



^ + (fc-l)A /i^ + Afc 



(3.7) 



fe-i 



3=1 



This permits us to conclude that 



H v +jX ^ + Xkj^\j- 



-IV'- 1 . 



3=1 



J 



V?-i. 



1 



ff + jX 



By inverting (3.8), we immediately arrive at the result (3.4). 
For v = 1, (3.8) can be written as 



(3.8) 

□ 



e-»y k (t)dt- 



e -At e -^t 



fc-1 

E 

3=0 



fc-1 

3 



(-l) j e- Xjt dt 



(3.9) 



e-^{e" A *(l-e- At ) fe - 1 }dt 

and this is an alternative derivation of the Yule-Furry linear birth process distribution. 

Remark 3.1. An alternative form of the distribution (3.4) can be derived by explic- 
itly writing the Mittag-Leffler function and conveniently manipulating the double sums 
obtained. We therefore have 



fc-i 



Pk 



fc-i 



Y(vm 



3=0 



fc-1 



(-At-)- v 

m=k K ' j=0 

(xtn k -\k-i)\ " 

T(v(k-1) + 1) 2->T(vm + l)^ 



fc-1 



fc-1 



(-iV(i + rr 
(-iV'(i + i) r 



(3.10) 



{-xry 



fc-i 



fc-i 
i 



(-iV'(i + i) r 
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The last step of (3.10) is justified by the following formulas (see 0.154(6) and 0.154(5) 
on page 4 of Gradshteyn and Ryzhik [6]): 

N , s 

^(-l) fe ( k ) (a + ky- 1 =0, valid for N>n>l, (3.11) 



k=0 

n 



k 

k I n 



k=0 



E^H fr (« + fc) n = (-!)"«!■ (3-12) 



What is remarkable about (3.12) is that the result is independent of a. This can be 
ascertained as follows: 

Sn=±(-D k (?) E ( H r ) = £,(") a'±(-lf (?) k-^-\ (3.13) 



r I * — ' \ r 

k=0 y ' r=0 y ' r=0 v ' k=Q 



By formula 0.154(3) on page 4 of Gradshteyn and Ryzhik [6], the inner sum in the third 
member of (3.13) equals zero for 1 < n — r +1 <n (that is, for 1 < r < n). Therefore (see 
formula 0.154(4) on page 4 of Gradshteyn and Ryzhik [6]), 



^ = (j)« E(-! 

^ ' k=a 



k I n 
k 



k n = (-l) n nl. (3.14) 



We now provide a direct proof that the distribution (3.4) sums to unity. This is based 
on combinatorial arguments and will subsequently be validated by resorting to the rep- 
resentation of N v (t) as a composition of the Yule-Furry model with the random time 

Theorem 3.2. The distribution (3.4) is such that 

oo oo k , , _ X 

£rf(*) = EE • - 1 i-^'^A-^n = i. (3.i5) 
fe=i fc=ij=i ^ ' 

Proof. We start by evaluating the Laplace transform L v (jx) of (3.15) as follows: 



m.)=ee 



fc=u=l 

oo fc— 1 



(3.16) 

k-l\ , „^ 1 



i — i „•_ n \ J / 
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A crucial role is played here by the well-known formula (see Kirschenhofer [8]) 

V ( N ) (-D k — = Nl (3 17) 

^ \ k J 1 ; x + k x{x + l)---(x + N)' 



k=0 



Therefore, 



/^"'v- (fc-1) 



^-i ~ r(; + i)iWA + i) 
a ^r( M VA + i + (/ + i)) 1 ' ' 



A 

i=o x " u i=0 

1 /-l roo 

(1 - aO"" 7 * -1 dx = / e^di, 
Jo 



A 

where B(p, q) = J Q a; p_1 (l — x) 9_1 dx for p, q > 0. This concludes the proof of (3.15). □ 

The presence of alternating sums in (3.4) imposes the check that p%(i) > for all k. 
This is the purpose of the next remark. 

Remark 3.2. In order to check the non-negativity of (3.4), we exploit the results of the 
proof of Theorem 3.2, suitably adapted. The expression 

E / e-^put)dt=^-j2 B [ k 'T +1 ) (3 - 19) 

k=l J ° fe=l ^ ' 

which emerges from (3.18) permits us to write 



oo fi ,,v—~y 

c 



dfc = I x^^-il-xf^dx 



1 x k-i^ e (y/x)Mi-) dx 
X 

A 

/ x k-l^ e -^Jl e -^ r /^r) dx _ 



(3.20) 
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The terms 

e -^x r /(Xr) = Ee -»Xr = I e -nt q lfat)dt (3.21) 



are the Laplace transforms of stable random variables X r = S(a r , 1,0), where a r = 



(v- cos z ^) 1 / !y (for details on this point, see Samorodnitsky and Taqqu [14], page 15). 



2A 

diffusion equation 



The term M 2A exp(— M J ) is the Laplace transform of the solution of the fractional 



~d^~ X d^ 1 o<i/<i, (3 22) 

u(x, 0) = S(x), 

with the additional condition that Ut(x, 0) = for 1/2 < v < 1, and can be written as 

u 2v {x,t) = — ± t^^ds (3.23) 

2AL(l-^)7 (i-s)" 

(see formula (3.5) in Orsingher and Beghin [10]), where Pv(x, 1) = ql(x, 1) is the stable 
law with <j\ = (j cos z ^) 1 / l/ . We can represent the product 



where 



/ e-" t q l/ {x,t)dt='[[c-^ xr ^ Xr l (3.25) 

•'O r=2 

Thus q^Xjt) appears as an infinite convolution of stable laws whose parameters depend 
on r and x. In the light of (3.24), we therefore have that 



/ e-»Y k (t)dt = 2 
Jo Jo 



»1 r t 



e -»t / x K-l / u , 2iy ( 2 ; jS ) giy ( a ; j t_ S )dsd.Tdt. (3.26) 







Since p^.(t) appears as the result of the integral of probability densities, we can conclude 
that pl(t) > for all k > 1 and t > 0. 

We provide an alternative proof of the non- negativity of p^(t), t>0, and of YlkPki^) ~ 
1, based on the representation of the fractional linear birth process N v (t) as 

N v (t)=N(T 2l/ (t)), <!/<!, (3.27) 
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where T2„(i) possesses a distribution coinciding with the folded solution of the fractional 
diffusion equation 

d 2v u d 2 u 

d^ = d^' o<^<i> (3 28) 

u(x, 0) = S(x), 

with the further condition that ut(x,0) = for 1/2 < u <1. 

Theorem 3.3. The probability generating function G v (u,t) = Ku Nv ^ of N v (t), t > 0, 
has the Laplace transform 

/ e-» t G u (u,t)dt= - —^"V *dt. (3.29) 

Jo l-u(l-e Af ) 

Proof. We evaluate the Laplace transform (3.29) as follows: 

/ e~^G v {u,t)At= e^V U fc V p" 1 (-rX^-AjOdi 
Jo Jo fc=1 J=1 \J V 

00 fe / 1 i \ i 

•:0(- l > 7 '~ 1 M 



fc=l 3=1 

7 J _ 1 OO fe— 1 



A*" + A? 



S^l*- 1 ^-!/—!— (M3.17)) 



A (^/A + l)(^/A + 2)..-0i'VA + fc) 



fe=i 



-T" ^ " (M"/A + l)-(/i"/A + l + ( } 

v-l /-l 00 



A 



/ ^w'a; i (l-x)^ /A dx (for < u.t < 1) 
;=o 



■dx=(l-x = e~ xt ) 



A Jq (1 — ux) 
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Remark 3.3. In order to extract from (3.29) the representation (3.27), we note that 

J ^ t ^u k Pv{N(T 2v (t)) = k}Ut 

{,.00 oo 
J ^2 ukpr { N ( s )^ k }fT 2 As,t)ds\dt (3.31) 

G(u,s)^ l/ - 1 e-^ s ds, 



which coincides with (3.29). It can be shown that 

/>oo 

/ e~ lit f T2 A^t)dt = n lJ - 1 e- sil \ s>0, (3.32) 
Jo 

is the Laplace transform of the folded solution to 

with the initial condition u(s,0) = <5(s) for < v < 1 and also u t (s,0) = for 1/2 < v < 1. 
In the light of (3.27), the non- negativity of p%(t) is immediate because 

/■OO 

Pi{N„(t) = k}= / Pr{iV(s) = fc}Pr{T 2l/ (f)eds}. (3.34) 
Jo 

The relation (3.34) immediately leads to the conclusion that X]fc°=i P r {-^W = &} = 1- 

Some explicit expressions for (3.34) can be given when the Pr{T2„(i) <E ds} can be 
worked out in detail. 

We know that for v = 1/2™, we have that 

Pr{T 1/2 „-i(t)ed*} 

= Pr{|Bi(|B 2 ("-|Bn(*)| •••)!)€ da} ( 3 - 35 ) 

/■oo e -s 2 /(4wi) /-oo e -wJ/(4w 2 ) /-oo e -w£_i/(4t) 

= ds2 ra / — do;i / — da;2 • • • / ; — da; n _i. 

J y'Anuji J \JAjxsjJ2 Jo V^nt 

For details concerning (3.35), see Theorem 2.2 of Orsingher and Bcghin [11], where the 
differences of the constants depend on the fact that the diffusion coefficient in equation 
(3.33) equals 1 instead of 2' 1 / 2 ) -2 . The distribution (3.35) represents the density of the 
folded (n — l)-times iterated Brownian motion and therefore B\, . . . ,B n are independent 
Brownian motions with volatility equal to 2. 
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For v = 1/3, the process (3.27) has the form Ny 3 (t) = N(\A(t)\), where A{t) is a 
process whose law is the solution of 



—-j, u(x,0) = 5(x). 



dt 2 / 3 dx 2 ' 

In Orsingher and Beghin [11], it is shown that the solution to (3.36) is 

3 1 



«2/a(aJ,<) = 2~yf t Ai \ 1757 )' 



'■it 



where 



Ai(x) 



cos ax 



da 



is the Airy function. Therefore, in this case, the distribution (3.34) has the form 



Pk W 



-As 



Ai 



ds, 



fc>l,i>0. 



(3.36) 



(3.37) 



(3.38) 



(3.39) 



Remark 3-4- From (3.3), it is straightforward to show that the probability generating 
function G„(u,t) = Eu N "^ satisfies the partial differential equation 

^G(u,t) = Xu(u-l)-^G(u,t), 0<i/<1, (3 40) 

G(u,0) = it, 

and thus EiV„(i) = ^| u=1 is the solution to 

^EiV„ = AE7V„ < i/ < 1, (3 41) 

EAT V (0) = 1. 

The solution of (3.41) is 

EN v (t)=E vA (Xt v ), t>0. (3.42) 
Clearly, the result (3.42) can be also derived by evaluating the Laplace transform 



-^EN^t) At 



J c P*{N(s) = k}Pi{T 2v (t)eds}\dt 

/ e""* / e Xs Pr{T 2 „(t)€ds}dt 
Jo Jo 

Jo n v — X Jo ' 
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and this verifies (3.42). The mean value (3.42) can be obtained in a third manner: 

L (-A^)c- A *di 



r°° . . / , \ /»oo 

/ e-^(i)^^ : (-ir 1 / E v 
Jo k=1 j=1 \J V Jo 



fe=i j=i 



A ^ (^/A + l)---(^/A + fc) 

^-i ~ r(fc)r^7A + i) 



(3.43) 



.y k i 

A T^/X + k + l) 

/ ^2kx k - 1 {l-xY"/ x ^-^ 7 —= e-^E Utl (Xt u )dt. 
J t, ^ • o 



A 



The result of Remark 3.4, EiV„(t) = i^i^At"), should be compared with the results of 
Uchaikin, Cahoy and Sibatov [16]. 

An interesting representation of (3.42) following from (3.27) gives that 

/•OO />OG 

EN v (t) = / e Xs Pi{T 2l/ (t)eds} = / EN (a) Pr{T 2 „(t) G ds}. (3.44) 
Jo Jo 

The expansion of the population subject to the law of the fractional birth process is 
increasingly rapid as the order of fractionality v decreases. This is shown in Figure 1 
and this behavior is due to the increasing structure of the gamma function for v > ap- 
pearing in the Mittag-Lcffler function E v ,\. This qualitative feature of the process being 
investigated here shows that it conveniently applies to explosively expanding populations. 



Remark 3.5. By twice deriving (3.40) with respect to u, we obtain the fractional equa- 
tion for the second-order factorial moment 

E{N v (t)(N v (t) - 1)} = g v (t), (3.45) 

that is, 

^g v (t) = 2Xg v (t) + 2XEN v (t), 0<i/<1, (g M) 
9,(0) = 0. 
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The Laplace transform of the solution to (3.46) is 



2/Li 



1 



1 



ji v - 2A n v — \ 
The inverse Laplace transform of (3.47) is 

E{iV„(t)(A^) - 1)} = 2£„, 1 (2Ai l/ ) - 2E^{\t v ). 
It is now straightforward to obtain the variance from (3.48), 

VarJV tf (t) = 2E„,i(2AO - ^,i(At") - ^(Af )■ 
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(3.47) 



(3.48) 



(3.49) 



For v = 1, we retrieve from (3.49) the well-known expression of the variance of the linear 
birth process 



VariVi(i) =e A *(e At - 1). 



(3.50) 



100 



50 



-■ 



V = 


1 


V = 


1/2 




1/3 




1/4 


v = 


l/sqrt(2) 



/ / 

/ / 

/ / 
/ / 



1 2 
Time (arb. units) 



Figure 1. Mean number of individuals at time t for various values of v. 
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Remark 3.6. If X\ , . . . , X n arc i.i.d. random variables with common distribution F(x) = 
Pr(X < x), then we can write the following probability: 

Pr{max(Xi,...,X JVi/(t) ) < x} 

CO 

= ^(Pr{X < x}) k Pr{N„(t) = fc} (by (3.27)) 



k=l 



= / G(F(x), S )Pr{T 2 „(t)eds} 
Jo 

f°° F(x)e~ Xs 

L i-F(;)ii-o->-) F ^'" £d "' 



(3.51) 



Analogously, we have that 

Pr{min(JCi,...,JCjv„(t)) > a;} 



(3.52) 



Remark 3.7. If the initial number of components of the population is no, then the 
p.g.f. becomes 

E(u N ^\N„(0) = n ) 

(3.53) 

= ^V+«o f°° e -Xzn ( n ° + ^~ 1 ^j (l- C - Xz ) k Pr{T 2l/ {t)e dz}. 
k=0 Jo \ / 

From (3.53), we can extract the distribution of the population size at time t as 

Pr{JV„(i) = k + n o \N v (0) = n } 



no + k — 1 
k 



(3.54) 

jj c- Xzna {l-c- Xz ) k Pr{T 2v (t) € dz}, fc>0. 



If we write fc + ?io = fc', then we can rewrite (3.54) as 
Pr{iV I/ (i) = fc / |iV !/ (0)=no} 



jfe'-l 
fc' — no 



(3.55) 

poo 

/ e - Az ™°(l-e- Az ) fc '~" o Pr{r 2t/ (0e dz}, fc'>n , 
Jo 



where fc' is the number of individuals in the population at time t. For no = 1, formulas 

(3.54) , (3.55) coincide with (3.4). The random time T2 V (t), t > 0, appearing in (3.54) and 

(3.55) has a distribution which is related to the fractional equation 

d 2u u d 2 u 

7^ = ^' 0<"<1- (3-56) 
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It is possible to slightly change the structure of formulas (3.54) and (3.55) by means 
of the transformation Xz = y so that the distribution of T2 V (t) becomes related to the 
equation 

|^ = A J g, 0<„<!, (3.57) 

where (3.1) shows the connection between the diffusion coefficient in (3.57) and the birth 
rate. 

Remark 3.8. If we assume that the initial number of individuals in the population 
is N v (0) — no, then we can generalize the result (3.4) offering a representation of the 
distribution of N„(t) alternative to (3.55). If we take the Laplace transform of (3.55), 
then we have that 

e -»t PT{N„(t) = k + n \N v (0) = n o }dt 

J™ + J\- Xzn °{l-e- Xz ) k Vr{T 2 „{t)^ dz}dt (by(3.32)) 

J™ ( no + ^ 1 ^je- Xzno (l-e' Xz ) k ^ 1 e-^ z dz 
"° + k~ 1 ) I" c- z{Xno+ ^\l ~ e~ Xz ) k dz 



o 



(3.58) 







n +k - 1 ^ J2 ^ J™ e - Z (\n +\r+»") dz 



k 



r=0 
k 



n + k—l\ „_jr^ [ k\ /_-.\r 
r=0 



r J \(n +r)+i^' 
By taking the inverse Laplace transform of (3.58), we have that 
Pr{N v (t) = k + n \N u (0)=n Q } 
rip + k — 1 \ v— v / k 



(3.59) 



r=0 



2(")(-l) r ^,i(-("o+r)An- 



From (3.59), we can infer the interesting information 
Pr{jV„(d<) = n + 1\N V (0) = n } 

n °E (l) (-lYE v> i(-(no + r)X(dty) (3.60) 



r=0 
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= n [E v>1 (-n X(dt) v ) - E v%l {-X{n + l)(dt)")] ~ n 



A(dt)" 

r(f + i) 



by writing only the lower order terms. This shows that the probability of a new offspring 
at the beginning of the process is proportional to ((It)" and to the initial number of 
progenitors. From our point of view, this is the most important qualitative feature of 
our results since it makes explicit the dependence on the order v of the fractional birth 
process. 

Theorem 3.4. The Laplace transform of the probability generating function G v (t,u) of 
the fractional linear birth process has the form 



H v (n,u) 



e-"*G„(t,ti) df = 



V-l fl 



X 



(1 -xY u i x 

1 — xu 



dx, 0<u<l, fi>0. (3.61) 



Proof. We saw above that the function G v solves the Cauchy problem 

0<v<l, 



d u G v . , dG v 
\u[u — 1) 



G„{u,0) = u. 



du 



(3.62) 



By taking the Laplace transform of (3.62), we have that 



// ' H v — /j," u = Xu(u — l)- 



du 



(3.63) 



By inserting (3.61) into (3.63) and performing some integrations by parts, we have that 
,2i/-i r 1 n _ r \n"/\ 



UpL 



X 



1 — xu 



■ dx — u/j, 



Xu(u — 1) 
: Xu(u — 1) 



u-i r 1 



1 — xu 



(i-xy u / x x 

A J (1 — xu) 2 



dx 



t-l A 



X 



■ dx ' 



A 



1 — xu 



1 — xu 



x=l 



,v-\ fl 



(1-x) 



A»" /A 



X J (1 — xu) 

uiu-i)/! 2 "- 1 f 1 x(i - xy/ x - 1 

A 



da; 



x=0 

- 1 xil-xY"/^ 1 



A 2 J (1 — xu) 



d.r 



(3.64) 



= — Uf-l 



2u-l /■! 



(1 — xu) 

(i - x y u / x 

X J Q (1 — xu) 
and this concludes the proof of Theorem 3.4. 



d.r 



dx, 



□ 
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Remark 3.9. We note that H u (fi,u)\ u= i = 1/fj, because G v (t,l) = 1. Furthermore, 



On 



,v-l 



\i v - A 



(3.65) 



which accords well with (3.42). 
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